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We use quantum process tomography to characterize a full universal set of all-microwave gates on 
two superconducting single-frequency single-junction transmon qubits. All extracted gate fidelities, 
including those for Clifford group generators, single-qubit tt/4 and tv/8 rotations, and a two-qubit 
controlled-NOT, exceed 95% (98%), without (with) accounting for state preparation and measure- 
ment errors. Furthermore, we introduce a process map representation in the Pauli basis which is 
visually efficient and informative. This high-fidelity gate set serves as another critical building block 
towards scalable architectures of superconducting qubits for error correction schemes. 



PACS numbers: 03.67.Ac, 42.50.Pq, 85.25.-j 

A critical prerequisite for building a scalable fault- 
tolerant quantum computer is the application of error 
correction codes [1-3]. In order to employ these codes, 
the underlying quantum gates must be performed with 
high fidelity above certain threshold levels. Experimen- 
tally, the onus currently lies in making scalable physical 
systems with universal quantum gate sets which surpass 
these thresholds. 

Two-dimensional error-correction surface codes, with 
gate fidelity thresholds of ^90-99.5% depending on mea- 
surement errors [2-4] , are particularly well-suited for su- 
perconducting qubit quantum processors, as repetitive 
tiling of qubit and resonator networks are apposite to 
proposed nearest-neighbor lattices [5]. Regarding su- 
perconducting qubits, single-qubit average gate fideli- 
ties exceeding 99% have been shown [6, 7], two-qubit 
gates are capable of entangled state fidelities greater 
than 90% [8, 9], and three-qubit entanglement has been 
observed [10]. Furthermore, recent experiments show- 
ing tenfold increases in coherence times of Joscphson- 
junction qubits [7, 11, 12] suggest that characterization 
of a complete high-fidelity universal set of quantum gates 
in superconducting qubits should be realizable. 

In this Letter we report gate fidelities greater than 
95% for a complete universal set of gates for two fixed- 
frequency superconducting qubits. The qubits arc single- 
junction transmons (SJT), coupled via a coplanar waveg- 
uide resonator [13]. The coherence times for our two SJT 
device are over twice as long as those of any previously 
reported superconducting multi-qubit system [9]. The 
gates characterized in this work include the single-qubit 
rotations {I, X T , X n/2 , Y n/2 , X w/i , l l/8 }® 2 (Re repre- 
sents a rotation of angle 9 around axis R) and a two- 
qubit CNOT. The combination of the 7r/4 rotation (com- 
monly referred to as the n/8 gate [14]), the Clifford group 
generators {^72,^/2}) an d the two-qubit CNOT form 
a universal set [1]. The two-qubit CNOT between the 
fixed-frequency qubits is implemented using the cross- 



resonance (CR) interaction [9, 15, 16]. Although the gate 
fidelities arc obtained via standard quantum process to- 
mography (QPT) of each respective gate, in this work we 
introduce a different, but efficient visual representation 
of process maps, the Pauli transfer matrix, 1Z. 7Z de- 
scribes the action of a process on the components of the 
density matrix represented in the basis of Pauli operators 
and helps establish a number of properties of the under- 
lying process which are otherwise hidden in the standard 
chi-matrix representation [14]. We find slight variations 
in the fidelities across the different gates and, from the 
associated maps, identify that our errors are not coher- 
ence limited but rather control limited and inherent in 
the QPT scheme. This suggests that as superconduct- 
ing qubits continue to improve in coherence times, the 
adoption of longer control sequences inspired by nuclear 
magnetic resonance quantum computing [17], and other 
efficient gate characterization schemes, will become nec- 
essary. 

The transmon qubit has become a popular choice for 
superconducting quantum computing applications due to 
its excellent coherence properties [13]. Implementations 
of transmons in planar qubit-cavity devices have tradi- 
tionally been formed via capacitively shunting a Cooper- 
pair box, preserving tunability of the qubit transition 
frequencies. The ability to dynamically tune qubit fre- 
quencies is necessary for a number of entangling gate 
schemes [8, 18] but often involves fast flux-biasing to fre- 
quencies with reduced coherence. In this work, we focus 
only on fixed-frequency SJTs, having previously gener- 
ated entangled states using an all- microwaves scheme [9] . 

Two SJTs [Fig. 1(a)] with transition frequencies 
wi/2tt = 5.0554 GHz and cj 2 /2tt = 4.9895 GHz are 
coupled via a 7.325 GHz coplanar waveguide resonator. 
Single-qubit rotations are performed by irradiating each 
qubit through independent on-chip microwave bias lines 
at the respective transition frequency. The relaxation 
times of the two qubits are measured to be = 8.2 fis 
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FIG. 1. (color online) Single-junction transmon and quan- 
tum process tomography pulse sequence, (a) Optical micro- 
graph of single-junction transmon device. Design is similar 
to standard transmon devices, with interdigitated capacitors 
shunting the Al/AlOx Josephson junction on either side. In- 
set: scanning-electron micrograph of Josephson junction, with 
critical current Iq = 28 nA. (b) Pulse sequence for performing 
quantum process tomography. Over-complete set of rotations 
{/, X„, -X±7r/2j ^±7172} are used to generate input states and 
to analyze the process before joint readout via the cavity. 

(2) 

and = 9.7 /_ts with Ramsey- fringe coherence times 
of T 2 (1) = 7.1 and T 2 (2) = 10.3 fjs. We attribute the 
improved coherence times to reduced surface loss contri- 
butions [19] via larger qubit feature sizes [10 fim capaci- 
tive shunt fingers and gaps, see Fig. 1(a)] and meticulous 
radiation shielding techniques [11, 20]. Further sample 
details are given in supplementary material [21]. 

We employ a series of repeated pulse experiments to 
accurately calibrate microwave amplitudes, offsets, and 
phases of single-qubit rotations [21]. The single-qubit 
rotations are shaped gaussian envelopes (gaussian stan- 
dard deviation o = 10 ns, total gate length 4cr) with 
derivative gaussian quadrature corrections to account for 
excited state leakage [6, 22]. 

The two-qubit interaction is the CR effect and its 
implementation and advantages as a multi-qubit gate 
scheme are previously detailed in Refs. 9 and 15. To 
review, the CR^ effect is exhibited as a qubit i (con- 
trol) state dependent drive of qubit j (target). In our 
system, we drive microwaves resonant with qubit 2 onto 
the microwave bias line addressing qubit 1. A residual 
classical cross-talk term is present, characterized by the 
parameter 777,12 = 0.22, which is the fraction of the direct 
Rabi frequency experienced by qubit 2 through driving 
on qubit 1. The strength of the desired two-qubit quan- 
tum effect is bounded by the residual qubit-qubit cavity- 
mediated dispersive interaction, which we determine to 
be J/ir — 7.44 MHz from detailed spectroscopy of trans- 
mon energy transitions and comparison to diagonaliza- 
tion of the Hamiltonian for our multi-level system. The 
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FIG. 2. Quantum process tomography for X^/4 <S> I repre- 
sented as the Pauli transfer matrix TZ. (a) Experimentally 
extracted TZ for X^/4 ® I with gate fidelity F g — 0.9687. 
To illustrate the action of TZ, an input state 1 00) (state fi- 
delity F s — 0.9818) is shown above TZ, and the output state 
cos(tt/8) J00) - isin(7r/8) |10) (F B = 0.9969) is shown to the 
right, (b) Ideal TZ for X w/i ® I. 



experimental tune-up and implementation of a CNOT 
gate with this interaction follows that of Ref . 9 and in this 
work the pulse-shape is a flat-top gaussian (a = 10 ns) 
of length 110 ns. 

To measure the two-qubit state, we employ a series of 
joint two-qubit measurements via non-linear driving of 
the cavity [23]. Calibrating this joint measurement [21], 
we perform state tomography on arbitrary two-qubit 
states through the application of an overcomplete set 
of analysis pulses {/, X^, X ±7T / 2 , Y ±7r / 2 }® 2 - An efficient 
maximum-likelihood estimation (MLE) algorithm [24] is 
used to compute two-qubit states which can be repre- 
sented by the Pauli state vector p [25]. The elements 
of p are the expectation values of the two qubit Pauli 
operators, (UV), where U,V £ I, X, Y, Z. 

QPT is accomplished through the compilation of an 
overcomplete set of 36x36 measurements [Fig. 1(b)]. 
State tomography as described above, is performed 
for the 36 different input states generated by 
{/, X v , X ±7I /2, Y ±7T / 2 }® 2 . Instead of the standard cm- 
matrix representation of the process map [14], we present 
the Pauli transfer matrix TZ, which maps an input Pauli 
state vector p- m to an output Pauli state vector p ut) 
Pout = TZp m . TZ is obtained through a semi-definite 
program taking into account the covariance matrix for 
the different independent measurements. Here, the semi- 
definite program extraction of TZ weights the measure- 
ments unequally [21]. The gate fidelity F g is calculated 
from the TZ map by F g = (Tr[Rlll] + d)/(d 2 + d) with 
d = 2n and n is the number of qubits. 

We perform QPT for the processes corresponding to 
the identity operation I® 2 , all independent single-qubit 
rotations of ir, ir/2 (around x- and y-axes), 7r/4, and ir/8, 
as well as the CNOT12 operation. The experimentally ex- 
tracted map TZ for the universal single-qubit ?t/4 rotation 
for qubit 1 is shown in Fig. 2(a) and the ideal 7?.idcai is 
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FIG. 3. Quantum process tomography for CNOT12 repre- 
sented as the Pauli transfer matrix TZ. (a) Experimentally 
extracted TZ for CNOT12 with gate fidelity F g = 0.9507. To il- 
lustrate the action of TZ, an input state (|0}+i | l))(|0}+i 1 1})/2 
(F s — 0.9787) is shown above TZ, and the output entangled 
state (| 00} + i |01> - |10) + i|H))/2 (F s = 0.9827, C = 0.994) 
is shown to the right, (b) Ideal TZ for CNOT12. 



shown in Fig. 2(b). From TZ, we estimate F g = 0.9687. 
The Pauli-transfer matrix TZ for the two-qubit entangling 
CNOT12 gate is shown in Fig. 3(a), from which we es- 
timate F g = 0.9507. Within this set the highest entan- 
gled state fidelity (F s — 0.9827) corresponds to the state 
(|00) +i|01) - |10)+j |ll))/2, with an associated concur- 
rence C — 0.994. The extracted F g for the complete set 
of gates are given in Table I. 

Figures 2(a) and 3(a) also demonstrate the action 
of the different Pauli transfer maps on specific input 
states. Above the experimentally extracted TZ are the 
measured Pauli state- vectors corresponding to 1 00) and 
(|0) + i |1»(|0) + i |l))/2 in Fig. 2(a) and Fig. 3(a), re- 
spectively. We can follow all the elements of each of these 
Pauli states downwards in the figures into TZ, which trans- 
fers these weights over into the Pauli operators to the 
right of TZ. Therefore, the Pauli state to the right re- 
flects the output state given the operation of the map 
TZ on the input Pauli state vector, and this visualization 
very simply demonstrates its effect. 

The TZ representation is a more efficient and intu- 
itive method for representing a quantum operation than 
the standard chi-matrix [14], as it consists of only real 
numbers and possesses a few other nice visual prop- 
erties. First, it is simple to tell if the map is trace 
preserving, which amounts to TZujk — &ij&ik f° r au 
j,k E {I,X,Y, Z}. Next, we can also determine if the 
map is unital, if TZjkji = SijSik for all j, k G {/, X, Y, Z}. 
Finally, the elements in TZ are bounded by ±1 and for any 
Clifford operation there is exactly one non-zero element 
in each row and column with unit magnitude. 

It is possible to further investigate the TZ maps to un- 
derstand the errors in our system and determine whether 
the loss in fidelity for all of the gates is due to statistical 
or systematic errors. For statistical errors, we use a boot- 
strapping method which generates new realizations of to- 
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TABLE I. Summary of gate fidelity F s , statistical error 
AF S , purified fidelity F pulc , and non-physical error T np = 
0.5||7?.mle — 72. e xp || 2 for complete universal set of gates on two 
qubits. 



mography experiments based off the variance of our mea- 
surement operator calibrations [21]. From these gener- 
ated experiments, new ensembles of maximum-likelihood 
estimates are obtained and the variance of these ensem- 
bles serve as an upper bound on statistical fluctuations 
of our estimated gate errors [21]. For all of the gates 
studied in this work, the statistical component is found 
to be ~ 3 — 6 x 10~ 4 (Table I), much smaller than the 
QPT extracted gate errors of s» 5%. 

The small value for statistical fluctuations suggests 
that our primary sources of error are systematic in na- 
ture. These can include decoherence, over- or under- 
rotations, and phase errors, all of which can occur dur- 
ing the actual processes to be characterized, or during 
the state preparation and measurement analysis gates. 
Given the coherence and gate times of the process to- 
mography sequence, we estimate a total error of 1.62% 
for single qubit gates and 2.55% for the CNOT, which is 
still smaller than the observed error and not the primary 
source. 

The remaining error can be accounted for via miscali- 
bration of 1% in the single-qubit gates used for prepara- 
tion and analysis. Simulations of the TZ maps including 
this level of calibration error agree well with the non- 
ideal elements found in the TZ data [21]. This is further 
confirmed through the calculation of the purified fidelity, 
-Fpurc, which estimates how close the unitary contribu- 
tion of the map is to the ideal gate. -F pU ro is defined as 
the overlap of the ideal map with the maximum eigenvec- 
tor of the Choi matrix [26] corresponding to 7\Lmle [21]. 
For all of the gates investigated, F purc > 0.99 (Table I) 
indicating that we perform the correct unitary to better 
than sa 1%, but we are adding on non-purity conserving 
operations due to systematic errors. These errors can be 
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Mechanism Error 

Statistical errors (AF S ) 0.1% 

Ti/Ta (simulation from 1 - F s ) 1.5 - 2.6% 

Unitary error in gate (1--F pure ) 0.5-1% 

Non physical errors ( T np ) 3% 

Total (1-F g ) 3.5-5% 



TABLE II. The error budget. 

the 2-norm distance between the raw experimentally de- 
termined lZ cxp and the most likely physical map 7?-mle- 
For all gates T np rs 3 — 4% (Table I), which suggests 
that propagation of small errors in our preparation and 
measurement operations in the QPT accounts for a large 
portion of the gate errors obtained experimentally [21]. 

Table II summarizes the different mechanisms to which 
we attribute our gate error budget. Further discussion re- 
garding these error metrics are given in the supplemen- 
tary material [21]. Our QPT results further reinforce the 
importance of other methods for gate characterization, 
such as randomized benchmarking [27] , as F s can be sig- 
nificantly misrepresented due to small systematic errors 
in preparation and analysis. 

Finally, as a precursor to full randomized benchmark- 
ing sequences, we show C and F B , after the repeated ap- 
plication of up to 12 CNOT gates on the input state 
generated by X_ 7T / 2 <8> /. Figure 4(a) reveals C for up to 
12 applications of CNOT, indicating as desired no entan- 
glement for even numbers (N) of gates. In Fig. 4(b), we 
plot the F s of the final state to the ideal states, the Bell 
state (|00) + i \ll))/y/2 for odd N and the input state 
(| 00) +i|10))/v / 2 for even N. By assuming an exponen- 
tial model for the state fidelity, F s = AF g + B, with A 
and B as fit parameters (dashed black line), we extract 
a gate fidelity F g = 0.9836. We find that this is in good 
agreement with a simulation of F g taking into account 
the coherence parameters of the system and the dura- 
tions of the gates involved in the experiment (solid blue 
line). 

Thus, we have shown a complete universal set of high- 
fidelity gates on two high-coherence fixed-frequency su- 
perconducting qubits. Although the gate fidelities ob- 
tained via QPT are at the 95% level, both the analy- 
sis of the 1Z maps with only 1% systematic calibration 
errors and the repeated CNOT sequences suggest that 
the intrinsic gate fidelities are > 98%, near the coher- 
ence time limit. Moving forward, finer pulse calibration 
tune-up sequences must be implemented and gate errors 
will be characterized via randomized benchmarking to 
avoid systematic errors in state preparation and analy- 
sis. Nonetheless, this demonstration further paves the 
road towards > 10 qubits for implementing more compli- 
cated pieces of quantum error correction. 

We acknowledge experimental contributions from J. 
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FIG. 4. (color online) Concurrence and state fidelity after N 
applications of CNOT. (a) Concurrence C of final states af- 
ter applying N CNOT gates to input state (|00> +i |10))/a/2. 
(b) Black squares are final state fidelity F s of experimentally 
obtained states to ideal Bell (original input) state for odd 
(even) TV. Dashed black line is a fit to the data assuming a 
model of AF^ 1 + B, where A, B, and gate fidelity F g are fit 
parameters. The fit gives an error per gate 1 — F s = 0.0164. 
Solid blue line is a theoretical model using the measured co- 
herence values and gate times. 
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SAMPLE PARAMETERS AND CALIBRATION 

The single-junction transmon ground to first excited- 
state frequencies are at f® 1 = 5.0554 GHz and f® 1 = 
4.9895 GHz. From spectroscopy of the system with vary- 
ing powers, we locate multi-photon transitions up to the 
third excited-state: / 1 02 /2 = 4.9454 GHz, f? 3 /3 = 4.8296 
GHz, / 2 02 /2 = 4.8753 GHz, = 4.7497 GHz. We es- 

timate the crosstalk parameter m\2 — 0.22 from Rabi 
driving qubit 2 via the bias line associated with qubit 1 
and comparing with the same drive but applied directly 
to qubit 2. From all the transitions of the two transmons, 
by using a model assuming transverse coupling between 
the two qubits, we estimate the coupling strength to be 
J/2tt = 7.44 MHz. We also estimate a residual ZZ inter- 
action of 1.1 MHz between the qubits via Ramsey fringe 
experiments with and without one of the qubits in the 
excited state. 

Single-qubit tt/2 and it pulses are calibrated us- 
ing sets of repeated sequences. First, the tt/2 am- 
plitudes are tuned up by measuring the cavity re- 
sponse after {X n /2>[X n /2,X n /2] N } for N up to 30. 
Then, we bootstrap off the tuned up n/2 pulse to 
find the appropriate ir amplitude via the pulse train 
{X„/ 2 , [^iD, again for N up to 30 pulses. Finally, 
the derivative-pulse shaping (DRAG) parameters [1] 
and gain imbalance between x and y rotations arc 
tuned via the experiments, {X v / 2 , [I T /2, —X n /2] N } and 
{X n/2 , [X n/ 2,Y w/2 , -Y w/2 , -X n/2 ] N } also for N up to 30. 

REPRESENTATIONS OF A MAP 

Defining Hd as a Hilbcrt space of dimensions d, and 
£(Hd) as the space of linear operators on Hd, then the 
Choi-Jamiolkowski isomorphism states that any map A : 
— > C(%d) can be represented by the Choi matrix 

PA = - d Y,E ij ®k{E ij ). (1) 

i-J 

Here pA G C(/Hd <8> 'Hd) and Eij is the matrix with 1 in 
the ij th entry and 0's elsewhere. In rcf. 2, Choi proved 
that A is completely positive if and only if P a is positive 
semidefinite. 

The Choi matrix allows us to test if the map is trace 
preserving by Tre[pA] = I/d and if the map is unital by 



TrA[/0A] = I/d. Furthermore, the action of the map on 
any state p G C(Hd) is given by 

A(p) =rfTr A [(p T ®/)p A ]. (2) 

An alternative representation of a map is the Pauli 
transfer matrix, 1Z, defined by 

TUj = ^Tr(P t A(P 3 )). (3) 

This representation is simply the Choi matrix expanded 
in the Pauli operator basis (here the Pj 's represent some 
ordering of the strings of Pauli operators {/, a x ,cr y , <r z }), 
since it is easy to show that 

Kij =Tr[p A ff (4) 

and therefore 

This is a useful representation as the act of the map on 
a state p in the Pauli representation (p = hp ■ P) is a 
simple matrix transformation of the Pauli state-vector p 
given by 

& = Kp, (6) 

which is obvious from our definition of 1Z. Therefore, the 
Pauli transfer matrix is also related to the Liouville or su- 
peroperator representation by a simple basis transforma- 
tion (the transformation from to P^) and composition 
of maps is described simply by matrix multiplication. 

The 1Z matrix has a few nice properties when compared 
to the standard chi-matrix [3] for determining whether a 
map is physical. First, 1Z consists of only real numbers 
and so contains exactly the same number of parameters 
as the gate, with no redundancy due to Hermiticity. Sec- 
ond, it is easy to see what operators are mapped to and 
from the identity by looking at the final row and col- 
umn respectively (assuming = I). If only the iden- 
tity is mapped back onto itself, then the map is trace- 
preserving, and if the identity remains the identity the 
map is unitial. This corresponds to a final row or col- 
umn of all zeros and a single one. Lastly, the elements 
of 1Z are bounded by ±1, and for any Clifford operations 
there is exactly one non-zero element in each row and 
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FIG. 1. A comparison of the theoretical and experimentally reconstructed 1Z matrices for the gates/ CgiX^/g, I®X n /4, I® X v /2, 
I <%> X n . 
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FIG. 2. A comparison of the theoretical and experimentally reconstructed IZ matrices for the gates X n / 8 ^>I, 1, 

X-x ® I. 
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FIG. 3. A comparison of the theoretical and experimentally reconstructed 1Z matrices for the gates I ® CNOT, ( 
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column with magnitude one, as Clifford operators sim- 
ply map Pauli operators to other Pauli operators. 
The average gate fidelity can be found by [4] 



_ dF p + l 
9 ~ d + 1 



(7) 



where the process fidelity is given by 

F p = Tr[p idcal p cxp ] = Tr[7^ cal -K cxp ]/d 2 , (8) 
assuming that the ideal map is a unitary map. 

QUANTUM PROCESS TOMOGRAPHY 
The Likelihood function 

Measurements in quantum mechanics can be rep- 
resented by positive operator valued measurements 
(POVMs). Mathematically, a POVM is a decomposition 
of the identity into positive scmidefinite Hermitian oper- 
ators X^fe=i Ek = I- The probability of each outcome k 
when state p is measured is Tr^E^p). In the continuous 
limit we have instead J dmE m = I and the probability 
density for measurement outcome m is p(m) = Tr(E m p). 

For a measurement of observable M — ^2 i m i \i)(i\ with 
Gaussian noise of variance v (as common with real am- 
plifiers) the POVM is 



-(ra-m,)7(2«) 



(9) 



(all units p.V). Having characterized the initial state and 
the measurement operator, we can completely determine 
1Z from data obtained for the applications of the pulses 
{/, X n , X ±7I /2, Y ±7T / 2 }® 2 before and after the operation 
to be analyzed. 

Defining m,j as the measurement result for the i th 
preparation and j th measurement we write the Likeli- 
hood function 



£(A)=n^iA)=n 



(12) 



where 



my = Tr[M,A(p 4 )] = ^Tr[p 4 P„]Tr[M,P m ]7e m „. 

ran 

(13) 

Here we have allowed the measurements to have different 
variances. Maximizing this likelihood function is equiva- 
lent to minimizing the negative log-likelihood function, 



Ao g = ^2,[m l3 - mij] 2 /vi 



(14) 



and so then the problem of Quantum Process Tomogra- 
phy (QPT) becomes 



minimize 8(A) T V 1 5(A) 
subject to pa > 0. 



(15) 



Here 5(A) = — and V is the covariance matrix. 



and the probability density for measurement outcome m 



is 



P(m\p) = -±= e- {m - mi)2/ ^P(i\p). (10) 



If we sample from this distribution N times and are inter- 
ested in fh — ^2 n mS n > /N then the central limit theorem 
asserts that as TV — > oo, the distribution of fh approxi- 
mates a normal distribution with mean Tr[Mp] and vari- 
ance v/N, 



P(fh\p) = 



1 



y / 2lTv/N 



-(m-Tr[Mp]) 2 /(2v/N) 



(11) 



To have complete information about the quantum op- 
erations we need to measure a set of linear independent 
observables and prepare linearly independent states. To 
achieve this we assume that the initial state can be reli- 
ably be prepared in the ground state and the measure- 
ment operator (M) is diagonal. Experimentally, from 
applications of {/, A^}® 2 pulses we reconstruct the out- 
comes of the measurement M. For this work we find 
M 00 = 0.00 35, M i = 0.0196, M w = 0.0302, and 
Mn = 0.0323 with a standard deviation of yjv = 0.0143 



Maximum Likelihood estimation using a 
semideflnite program 

A general scmidefinite program (SDP) can be ex- 
pressed in the standard form: 



minimize b x 



subject to Fq + FiXi>0 



(16) 



where the n-element real vector x contains the variables 
to be optimized with respect to the linear objective func- 
tion described by an n-element real vector b and con- 
strained by the constraints are described by a set of m- 
by-m Hermitian matrices F , Fi where i = l,...n. 

Since the objective function of Eq. (15) is quadratic, it 
is not in the standard form. In order to rectify this we in- 
troduce a slack variable t defined by t — 5(A) T V~ 1 8(A) > 
0. By the Schur complement the previous inequality is 
satisfied if and only if 




> 



(17) 
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TABLE I. A compilation of various error measure and measures of physicality for our universal set of quantum gates. The 
columns, from left to right, contain: the name of the gate, the process fidelity, the gate fidelity, the error in gate fidelity found 
from bootstrapping, the value of the Pauli transfer matrix from the identity to the identity (should be unity for a physical map), 
the purified fidelity, the maximum eigenvalue of the reconstructed state, the sum of the negative eigenvalues of the reconstructed 
state obtained by simple matrix inversion and finally three two norm distances between the pairwise combinations of the ideal 
gate (pideai), the physical maximum likelihood estimate (pmle) and the unphysical gate obtained by inverting the data (pdata)- 



and V > 0. Therefore, we can express the QPT problem 

as 



-Hi 

minimize b y 

subject to A = | | > 



(18) 



where b = {1,0^} and y = {t,r} with r being the 
column-major vector representation of 1Z. Although 
there are many numerical libraries for solving semidef- 
inite programs in this standard form, we choose to use 
SeDuMi [5]. 

In Fig. 1 - 3 we show the experimental 1Z ma- 
trix after applying MLE as well as the theoretical 
1Z. These are plotted for the measured operations 
{I,X n/8 ,X w/4 ,X n/2 ,X w ,Y n/2 }® 2 and CNOT. The cor- 
responding gate fidelities are shown in Table I, and are 
all above 0.95. 



ERROR ESTIMATION 
Statistical analysis of errors 

The conceptually simplest way to obtain a confidence 
region on our QPT extracted error rates would be to 
calculate 



P 



Is dpA£(p A ) 



dp A -P(p A |data) - f rf . 
5 Jv QPaApa) 



(19) 



where we interpret S as the region that contains the ac- 
tual state with probability p (V is the entire space of 



Choi matrices). Even if we model the noise on our mea- 
surements as Gaussian, constraining the domain of the 
likelihood function to the set of positive semidefinite ma- 
trices makes this integration extremely difficult. 

The simplest way to bound the error due to any noise 
process would be to take many independent realizations 
of an experiment. However, given reduced measurement 
fidelities, it can be costly to run the experiment a large 
number of times, and instead we can use a bootstrapping 
method to approximate a new realization of the noise on 
an existing data set. Given the measurement record fhij 
we can simulate 



- (fc) 



mi 



+ sJv j /NW^\ (20) 



(k) 

where is a gaussian random variable of mean and 
variance 1 and N — 10000. From this new measurement 
record we can compute an ensemble of maximum like- 

(k) 

lihood estimates Pmle- This ensemble is not the same 
as averaging over different realizations (in particular the 
mean is offset) but the variance of the distribution serves 
as an upper bound on the statistical fluctuations in the 
estimate. In Table I the error for one standard deviation 
gives about 5 x 10~ 4 , which is much smaller than the 
measured error (w 5%). 



Systematic errors 

There are a number of errors that can lead to sys- 
tematic errors such as decoherence, over-under rotations, 
phase errors, etc. These errors can occur in both the gate 
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FIG. 4. The TZ matrix for (a) a X n / 2 ® I with a 1% under rotation and 1% phase error, (b) I ® with a 1% under rotation 
and 1% phase error, (c) a amplitude damping channel (with error probability of 0.5), (d) a dephasing channel (with error 
probability of 0.5). 



being characterized as well as the gates used in the anal- 
ysis and input state preparations. Some examples of TZ 
maps for these types of errors are shown in Fig. 4. 

In our experiment the coherence times were = 
8.2 fjs, T[ 2) = 9.7 /US, T 2 (1) = 7.1 ^ts and T 2 (2) = 10.3 ^s 
and we can estimate that this would lead to an error of 
1.62% for the single qubit gates (total time for prepara- 
tion, process and analysis pulses of 120 ns) and 2.55% 
for the CNOT (total time 190 ns). This is much smaller 
than the observed error, making it likely that there are 
other, larger sources of error as well. This is confirmed by 
looking more closely at the TZ matrix, Fig. 1 and Fig. 2. 
The features associated with amplitude damping are not 
present, whereas we clearly see signatures of under rota- 
tions and phase errors, see the TZ plots of the X n / 2 ® / 
and / <E) in Fig 4. We estimate that this coherent 

error accounts for about 1% of the total gate error in 
our experiment. This is confirmed to some extent by the 
high values in the purified fidelity The purified fidelity 
is the overlap of the largest weighted eigenvector with 
the ideal state (-F pU re = Tr[pidealPA mal ]- This confirms 
that we are making the correct unitary with rj 0.5 — 1% 



error but we are adding on to it some non-purity con- 
serving operations from systematic errors (the weight of 
the maximum eigenvalue is about 0.94-0.96). 

To have a measure of the systematic errors we recon- 
struct a TZ matrix without enforcing the Choi matrix to 
be positive semidefinite. This is done by using 

f = {W T W)- 1 W T m (21) 

where rh is the column-major vector representation of m 
(measured data) and W is the transfer matrix defined 
by Eq. (13). To quantify how much the TZ matrix has 
changed due to the SDP, we calculate the total weight of 
negative eigenvalues (see Table I). Here we see that about 
10% of the eigenvalue weight is negative. Since rotational 
and phase errors in the operation pulse can not result in 
non-physical operations and our statistical fluctuations 
are small this must arise from errors in the preparation 
and mesaurement pulses. To estimate the error induced 
by this we calculate the 2-norm distance, 

\\P1-P2h = V^l\Pi-P2\] (22) 

between the Choi matrices of the maximum likelihood es- 
timate, a direct (nonphysical) inversion of the data, and 
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TABLE II. Simulations of tomography results when faulty 
preparation and measurement pulses are used to reconstruct 
a perfect identity operation. Each simulation is averaged over 
ten random gate set errors. 

the ideal map. These results suggest that we have about 
a 5% error from the ideal operation and the physical oper- 
ations have about a 3 — 4% error from the reconstructed. 
This is all consistent with the analyzing pulses having 
over/under rotations and phase errors which can be very 
loosely thought of as a mechanism that reduces the pu- 
rity of the operation (add more non-zero eigenvalues to 
the Choi matrix). 

One simple test of this conjecture is to simu- 
late the tomography data for a perfect process us- 
ing preparation and measurement pulses drawn from 
a faulty gate set. We take the 36 element gate set 
{/, Xn, X ±7r / 2 , Y±tt/2}® 2 an d apply an independent uni- 
tary error U e = cxp(— ieiJ ran d/2) to each gate, where 
-ffrand is a random Hamiltonian of unit norm. The re- 
sulting error in the gate fidelity for each preparation or 
measurement pulse is approximatly (e/2) 2 . In table II 
we simulate measurement data assusming a perfect iden- 
tity gate. Small errors in the analyzing gates can lead 
to significant errors in the tomgraphy. In particular, er- 
rors of around 1% lead to the expected 3 — 4% errors in 
reconstruction, even if the gate being analyzed is a per- 
fect identity transformation. This error due to process 
tomography appears to become more severe as a ratio 
of the single gate error as the gates get more precise, 
and we conjecture that it will also scale poorly with the 
dimension of the system. 
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